Tunable joint measurements in the dispersive regime of cavity QED 
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Joint measurements of multiple qubits open new possibilities for quantum information processing. 
Here, we present an approach based on homodyne detection to realize such measurements in the 
dispersive regime of cavity /circuit QED. By changing details of the measurement, the readout can 
be tuned from extracting only single-qubit to only multi-qubit properties. We obtain a reduced 
stochastic master equation describing this measurement and its effect on the qubits. As an example, 
we present results showing parity measurements of two qubits. In this situation, measurement of an 
initially unentangled state can yield with near unit probability a state of significant concurrence. 

PACS numbers: 03.65.Yz, 42.50.Pq, 42.50.Lc, 74.50.+r 



In most of the current quantum information exper- 
iments, measurements are used to extract information 
only about single-qubit properties. Joint measurements 
where information about both single and multi-qubit 
properties can be obtained offer new possibilities. Ex- 
amples are the test of quantum paradoxes [l[ , test of 
quantum contextuality [2j], realization of quantum state 
tomography with weak measurements [3J45J and cluster- 
state preparation A particularly powerful type of 
joint measurement is parity measurement, where infor- 
mation is gained only about the overall parity of the 
multi-qubit state, without any single-qubit information. 
This type of measurement can be used for the genera- 
tion of entanglement without unitary dynamics 
for quantum error correction [ill . and determinis- 
tic quantum computation with fermions fl3L \\a\ . In this 
paper, we show how such joint measurements can be real- 
ized in the dispersive regime of cavity QED (l5| . In par- 
ticular, we show how the character of the measurement 
can be tuned from purely single-qubit to parity read- 
out. As a realistic example, we present results for circuit 
QED [3- 5] and show that states with large concurrence 
can be obtained. Entanglement generation by measure- 
ment was previously studied in this system IB - HI, but 
ignoring information about the parity. With parity mea- 
surements, entanglement generation by measurement can 
be deterministic rather than probabilistic. 

We consider a pair of two-level systems (i.e. qubits) of 
frequencies uj aj with j = 1,2 coupled to a high-Q cavity 
of frequency co r . In the dispersive limit, where |Aj| = 
\(u aj — <*v)| >• \gj\ with gj the coupling strength of qubit 
j to the cavity, the Hamiltonian of this system takes the 
form [H 



This result is valid to second order in the small param- 
eter Xj — gj/Aj. Here, we have defined the dispersive 
coupling strength x 3 = gj^ji the Lamb-shifted qubit fre- 
quency ib aj and the strength of qubit-qubit coupling me- 
diated by virtual photons J q = (72(1/^1 + 1/^2)/2 [2(il |. 
The last term represents a coherent drive on the cavity 
of amplitude e m (t) and frequency u m f=s u> r , appropriate 
for measurement of the qubits. With this choice of drive 
frequency, we have safely dropped a qubit driving term 
of amplitude \jt m (20| . In order to focus on entangle- 
ment generated by measurement only, we drop the term 
proportional to J q . This is reasonable since the possi- 
ble measurement outcomes are eigenstate of the flip-flop 
interaction o-\a\ + o 1 _o\, as will be clear below. 

Coupling to unwanted degrees of freedom is mod- 
eled by using a Lindblad-type master equation (2l| . In 
Ref. [161 ]. a master equation for the qubits only was ob- 
tained by enslaving the cavity to the qubit dynamics. 
This approach is valid only in the limit where damping 
of the cavity k greatly overwhelms the dispersive cou- 
pling strength x 3 ■ Here, we go beyond these results by 
using a polaron-type transformation to trace-out the cav- 
ity 123. [2^ . Starting from Eq. fT}, we find following 
Ref. [22j the effective master equation 



(2) 

where T>[c\- = c • c) — {c*c, -}/2. In this expression, 7^ is 
the relaxation rate of qubit j and 7^ its pure dephasing 
rate. The fourth term represents Purcell damping at the 



while the last contains both measurement- 



+ e m (t)(a t e- l "'" t + h.c.). 



(1) 



rate \jn 

induced dephasing (I^ y ) and ac-Stark shift (A^ v ) by the 
measurement photons [? ]. In Eq. (J2|), x (y) stands 
for one of the four logical states ij with i,j <G {g,e} 
the qubit 's ground and excited states and = |rc) (rc[. 
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FIG. 1. (Color online) a) Phase space illustration of the sta- 
tionary states \otij) for: gi — —g-2 = — 15ac and x 3 ~ 1.5k. 
The drive is at resonance with the bare cavity A r = and its 
amplitude is e = k/2. b) Normalized rates of information gain 
for <f> = or = tt/2: T 01 {<f>) = T 10 (cf>) (full red line), r n (0) 
(dashed purple line). The vertical lines are at ±2x J ■ c) Nor- 
malized rates at A r = as a function of x 1 '■ Toi(0) = rio(0) 
(full red line), rn(7r/2) (dashed purple line). The vertical 
line indicates the value of x 1 used in panel b). Other rates 
Fn(O), roi(7r/2), and rio(7r/2) are zero and not shown. 



Measurement-induced dephasing and ac-Stark shift are 
given by 



T7 = {Xx - Xy)^[axa* y ], 



A x c v = (Xx -XyWa x a* y ], 



(3) 



(4) 



where Xx = (x\ X^- x J(J l\ x } an d ot x the amplitude of the 
coherent state when the qubits are in state \x). This 
amplitude satisfies 



-i{ui r 



Xxja x 



Ka x /2. (5) 



The reduced master equation Eq. ^ is a very good ap- 
proximation to the full dynamics when k/2> 7y. Since 
7ij does not include Purcell damping, this inequality is 
easily satisfied with current Purcell limited qubits [24|. 

To go beyond information about average evolution, we 
use quantum trajectory theory of homodyne measure- 
ment on the transmitted cavity field to obtain informa- 
tion about single experimental runs (25| . Following the 
approach of Ref. [22j , we find in the multi-qubit case the 
reduced stochastic master equation (SME) 

pj =Cpj + M[c4,]pjt(jt) - i[c+__ n/2 , pj]f(t)/2, (6) 

and the measured homodyne current is proportional to 
J{t) = Tr[c p, 7 ]+£(i). Here A4[c]- = {c, -}/2-Tr[c-]-, and 
£(t) is Gaussian white noise satisfying = and 

E[£(t)£(t')] = 6(t - t'), with E[] denoting an ensemble 



average over realizations of the noise. This stochastic 
equation is valid for k/2 3> 711 + 712, which is again 
easily satisfied [24j|. 

In Eq. the joint measurement operator is 



c = vT lo (0)oi + y/Toi(4>)o* + y/TuMalal , (7) 



where 

Tij((j)) = K,r)\j3i. 



(3ij = {a ee + (-l) j a eg + (-lYa ge + {-l) l+ >a gg ) /2, 



(8) 



with </> the phase of the local oscillator, 8 a = Arg(a), and 
77 the efficiency with which the photons leaking out of the 
cavity are detected, represents the rate of informa- 
tion gained about the first qubit polarization (zj=10), 
second qubit polarization (ij=01) or the parity (ij = 11). 
An optimal measurement occurs when 0^-^/2 = since, 
in this case, all the back-action arising from the measure- 
ment is associated with information gain (22j. Given the 
form of c ( p_ 7T /2, this cannot be realized, except in trivial 
cases. 

Given that x J > — — w m and <f> can be changed 
in-situ 043' the form of the measurement operator 
can be tuned (in the dispersive approximation, changing 
e m only leads to an overall rescaling). There are several 
useful choices of c^. For example, an equally weighted 
joint measurement (all |r,j| equal) is ideal for quantum 
state tomography since in this case both the required 
single and two-qubit information are on an equal footing. 
In the limit Ix 1 ± \ 2 \ ^ K , this is achieved by choosing 
uj m to match one of the four pulled cavity frequencies 
w r + Xx- As can be seen in Fig. [T})), for x 1 = X 2 1 an 
equally weighted joint measurement is realized by setting 
A.,. = ±2x' J . For this choice of x J however, at A r = it 
is not possible to determine which qubit is excited and as 
a result the measurement is either completely collective 
(a\ + a\) for <ft — [l(| or more interestingly extracts 
information only about the parity [a\a 2 z ) of the combined 
two-qubit state for = 7r/2. 

This can be understood by considering the steady-state 
cavity amplitude a x . Fig. |TJl) shows a phase-space plot 
corresponding to the four coherent states \a x ) for the pa- 
rameters given in the caption. Since x 1 = X 2 > the coher- 
overlap, while Im[a ee ] = Im[a gg ] 



ent states a eg and a 



but Re[a ee ] 7^ Re[a ffs ] [? ]. As a result, measurement 
of the Q {<j) = 7t/2) quadrature reveals information only 
about the parity and / (</> = 0) the collective polarization. 
Since, for these parameters, there is information in the 
quadrature orthogonal to the measurement, C0_ 7r / 2 7^ 0, 
and this measurement is not optimal. As illustrated in 
Fig. [Tp) however, as the ratio x J / K i s increased, the 
measurement becomes optimal for parity with Foi(0), 



normalized by £y=oi 10 
MX 3 ) 2 - 



n [r i: ,-(0) + Tij(ir/2)], scaling as 
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FIG. 2. (Color online) Mean fidelity to \<j>+) (dashed green 
and full purple lines) and |?/> + ) (dotted-dashed red and dot- 
ted blue lines) obtained from solving Eq. |[2j. Dashed green 
and dotted-dashed red lines: (jij, 7</>j) — (0,0). Full purple 
and dotted blue lines: (k/250,0). The measurement drive is 
etanh(i/(j) with a = 1/k and the initial state (|<?) + |e)) ® 
(\g) + |e))/2. The other parameters are g\ — —g-z — —100k, 
X j = 10k, e = k and k/27t = 5 MHz. 



An application of parity measurements is the genera- 
tion of entangled states from separable ones I7h10I. In 
contrast to collective polarization measurements [lq - ll ! , 
this can be achieved with unit probability. For example, 
with the initial separable state (\g) + |e)) ® (\g) + \e))/2, 
the measurement ideally projects on the Bell states 
10+) = (\eg) + \ge))/V2 or |^+) = (\gg) + \ee))/V2. That 
is, evolution under Eq. ([5]) shows a collapse of the separa- 
ble state to \4>+) or \ip+), conditioned on the record J{t) 
being predominately negative or positive respectively. 

There are four main causes of errors in this collapse. 
The first is relaxation and damping [the dissipative terms 
of Eq. ((2)]. Interestingly, with the parameters of Fig. [TJ 
Ai = — A2 such that \cj>+) is immune from Purcell de- 
cay [l(| . The second is the time-dependent ac-Stark shift 
[unitary contribution from the last term of Eq. @] which 
causes a phase accumulation between \gg) and |ee) in 
This contribution can be seen as a slow oscilla- 
tion of the fidelity F = (ip\p\ip) between the state 
and those obtained by numerical integration of Eq. ([2} . 
This is illustrated in Fig. [5] There, the mean fidelity to 
\4>+) is always 1/2, since half the density matrices col- 
lapse to that state, while oscillations due to the ac-Stark 
shift appears in the fidelity to However, this shift 

is deterministic and can thus be undone. The third er- 
ror comes from Co ^ causing a stochastic phase be- 
tween \gg) and \ee) [last term of Eq. ©]• For a given 
experimental run, this does not reduce the concurrence 
or purity of the state [because is known from J(t)]. 
However, since this phase varies from shot to shot, the 
ensemble averaged state is mixed. This error can be over- 
comed by performing J(rj)-dependent single qubit phase 
operations after the measurement or, more simply, by op- 
erating in the large x° limit where its effect is negligible 
as illustrated in Fig. [TJ;) . Finally, the measurement is not 
ideal in the sense that measurement-induced dephasing 
affects the measurement outcome \ip+) (i.e. T e ^' 99 7^ 0). 
However, this effect can be made negligible by increasing 
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FIG. 3. (Color online) a) Concurrence as a function of time 
for 10 4 trajectories. Dotted-dashed red line: (7ij,70j,»?) = 
(0,0,1), dashed green line: (k/250,0,0.8), dotted blue line: 
(k/250,0,0.2) and full purple line: (k/250,0,0.05). <f> = tt/2. 
All other parameters are the same as in Fig. [2] b-c) His- 
tograms of the integrated current s{i) at b) t = 1.6/k and c) 
t = 6.3/k. 



the ratio X j /n since rii(7r/2)/T"' 99 - {x* / n) 2 . 

To show that the system collapses to \4>+) or \(j>+), 
and that entanglement is generated with unit probability, 
Fig-El 1 ) shows the mean concurrence E[C(pj)], averaged 
over 10 4 trajectories. There, departure from unit con- 
currence in the dotted-dashed red line (no damping, unit 
detector efficiency 77 = 1) is only due to measurement- 
induced dephasing. The dashed green, dotted blue and 
full purple lines take into account relaxation with tt/^yij ~ 
250, and detection efficiency of 77 = 4/5, 77 = 1/5,77 = 
1/20, respectively. The latter corresponds to current ex- 
perimental values [24|]. The ratio n/jij is slightly out of 
reach of current experiments when taking into account 
that x 3 = 10k is also required. This cannot be achieved 
with transmons as current experiments have reached the 
maximal possible coupling [26j. However, new ideas to 
increase the qubit-cavity coupling can help in achieving 
these parameters [27]|. 

Small detection efficiency reduces the ratio Tu/T^' 99 , 
which in turns corrupts \ip+). As illustrated in Fig. [3ja,), 
this results in lower concurrences when 77 < 1. Interest- 
ingly, \<t> + ) is not affected by this detection efficiency [lit . 
Nevertheless, improvement in detection efficiency is re- 
quired to match concurrences that can be realized with an 
entangling Hamiltonian Recent improvements with 
near quantum-limited amplifiers are a good step in this 
direction (2j j. 

Having generated one of the two orthogonal entangled 
state, it is necessary to distinguish them efficiently. Us- 
ing the experimental record J(t) to compute pj(t) from 
the SME Eq. © is not efficient since the record is widely 
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FIG. 4. (Color online) Average concurrence (dashed red line), 
average fidelity to \<j>+) and \tp+) (dotted green line) and suc- 
cess probability (full blue line) as a function of the threshold 
at time t = 18.5/k. (7ij,7«,»?) = (re/250,0,0.05). The other 
parameters are the same as in Fig. [2] The ac-Stark shift af- 
fecting \ip+) has been corrected before evaluating the fidelity. 



In conclusion, we have shown how measurements in 
the dispersive regime of two-qubit cavity QED can be 
tuned from accessing single to multi-qubit information, 
thus allowing for example parity measurements. In ad- 
dition to allowing complete characterization of the two- 
qubit states 043 and the implementation of quantum in- 
formation protocols 0, [ll], , this allows for generation 
of entanglement by measurement with unit probability. 
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fluctuating. As a result, a useful and more efficient quan- 
tity to distinguish the states is the integrated current 




where rjj is the steady-state value of r n (7r/2). Fig.JSJs) 
and c) show two histograms of s(t) at times t = 1.6/k 
and t = 6.3/ k. These results are for r/ = 1 and exclude 
damping for illustration purposes. The full blue lines 
are Gaussians fits to the histograms. These separate at 
a rate ~ Tn. At times large compared to 1/Tn, but 
short compared to T\ and T2, the distributions are well 
separated and correspond to \ip+) and | </>+). 

As shown in Fig. [3t), we introduce a threshold sth to 
distinguish these states. All outcomes with s(t) < sq — s^ 
(condition c = — ) are assigned to | </>+), while those with 
s(t) > so + sth (condition c = +) to \ip+) where sq is 
the median of s. Values outside this range are disre- 
garded. A success probability P s can then be defined as 
the probability for s to be outside the range sq ± Sth- To 
quantify the success in generation and distinguishabil- 
ity of the entangled states, we define the average fidelity 
F= [{4> + \E-[pj]\4> + ) + (iP + \E + [pj}\tP + }]/2 and average 
concurrence C = [C(E+[pj})+C(E-[pj})}/2. E c [pj] rep- 
resents the ensemble average over pj for condition c = ±. 
These quantities are illustrated as a function of s t h for 
the fixed integration time t = 18.5/k in Fig. [4] Even 
when keeping all events (s t h = 0), F and C are large 
with values 0.92 and 0.79 respectively. That is, with this 
procedure, it is possible to create and distinguish highly 
entangled states with unit probability. If willing to sac- 
rifice some events, this average fidelity and concurrence 
is increased to 0.98 and 0.91, respectively. The deviation 
from unity in the large s t h limit is due to slight corruption 
of the state discussed previously. 
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